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Abstract 

We address the problem of a non-perturbative formulation of superstring theory by means 
of the recently proposed matrix models. For the model by Ishibashi, Kawai, Kitazawa and 
Tsuchiya (IKKT), we perform one- loop calculation of the interaction between operator- like 
solutions identified with D-brane configurations of the type IIB superstring (in particular, for 
parallel moving and rotated static p-branes). Comparing to the superstring calculations, we 
show that the matrix model reproduces the superstring results only at large distances or small 
velocities, corresponding to keeping only the lowest mass closed string modes. We propose a 
modification of the IKKT matrix model introducing an integration over an additional Hcrmitian 
matrix required to have positive definite eigenvalues, which is similar to the square root of the 
metric in the continuum Schild formulation of IIB superstrings. We show that for this new 
matrix action the Nambu-Goto version of the Green-Schwarz action is reproduced even at 
quantum level. 
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1 Introduction 



Recent interest in matrix-model formulation of superstrings has been initiated by the proposal 
of Banks, Fischler, Shenker and Susskind [jlj that the non-perturbative dynamics of M theory is 
described by supersymmetric nx n matrix quantum mechanics in the limit of large n. This matrix 
model has been investigated in a number of subsequent papers pi |4|, ||, n], and its operator-like 
classical solutions were identified with D(irichlet) p-branes (for even p) of the type IIA superstring 
theory. 

Another matrix model has been proposed by Ishibashi, Kawai, Kitazawa and Tsuchiya 
(IKKT) for type IIB superstrings in ten dimensions. The action of the model is defined by 

S = a ^Tr[A fl ,A u ] 2 -^Tr(^[A^^ + f3n , (1.1) 

where A^ and ip a are n x n (n — > oo) Hermitian bosonic and fermionic matrices, respectively. 
The parameter n is considered as a dynamical variable which makes a crucial difference between 
the action (|l]l]) and the one of (dimensionally reduced) ten-dimensional super Yang-Mills. The 
action ( |l.l| ) is associated with the IIB superstring in the Schild formalism (with fixed K-symmetry): 

Sgchiid = / d 2 a L(^{X* t X v } 2 PB - % -^{X^} PB )+f3^ , (1.2) 

where the commutators are substituted by the Poisson brackets. Properties of the IKKT matrix 
model were further studied in [||, ^, 10|.El 

The Dp-branes with odd p of type IIB superstring theory appear in the matrix model ( |1.1| ) as 
operator-like solutions of the classical equations 

[A", LV A v ]] = , [Ap , (r*V)«] = . (1.3) 

A general multi-brane solution has a block-diagonal form and is built out of single p-branes. The 
solution associated with one p-brane is given by 

Aj=(Pi,Qi,...,P e} x,Q 2 +i,0,...,0) J C 1 = 0, (1.4) 

where P's and Q's form (p+ l)/2 pairs of operators (infinite matrices) obeying canonical commuta- 
tion relation on a torus associated with compactification (of large enough radii L a /27r) of the axes 
0, . . . ,p so that L a L a+ i/n 2 /( p+1 ) (even a) is kept fixed as n — > oo. This solution for D-string (p = 1) 
was constructed in (7) in analogy with ||] and was extended for p > 3 in |9|, [l0| in analogy with (5| . 
One of the arguments in favor of this construction is based on the correct large-distance behavior 
of the one-loop matrix model calculation of the interaction between anti-parallel D-strings @ and 
higher p-branes |J. 

In the present paper we continue the investigation of the IKKT matrix model. We perform 
comprehensive one-loop calculation of the interaction between two Dp-branes and compare the 
matrix-model and superstring results. For parallel moving p-branes we demonstrate that the matrix 



model reproduces Bachas' superstring result [13] only at large distances or small velocities keeping 



only the lowest mass closed string modes. We propose a modification of the IKKT matrix model 



1 Another approach to the type IIB superstring is discussed in O, ^2). 
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introducing an additional (positive definite) Hermitian matrix Y %3 , which is a dynamical variable 
to be identified with ^J~g in eq. (|1.2| ), The classical action has the form 

S d = a (-1 TrF- 1 ^,^] 2 - ~ Tr(^L4^])) + /?TrY (1.5) 
and reduces to the classical action of the non-abelian Born-Infeld (NBI) type: 



a 



51 = y/afi Tr ^-{A^A v f - - Tr (^[A^, rj,]) , (1.6) 

using classical equation of motion for Y. Moreover, we show that the Nambu-Goto version of 
the Green-Schwarz action is reproduced even at the quantum level, if one chooses appropriately 
the measure of integration over Y. Our results are therefore much more general than displayed in 



eqs. (|E5D and flOp . 

In Sect. 2 we perform one-loop calculation of scattering of parallel p-branes in the IKKT matrix 
model (as well as the interaction of rotated static p-branes) and compare to the superstring calcu- 
lations. We also reproduce previously known results for the anti-parallel p-branes using the new 
technique. In Sect. 3 we propose a modification of the IKKT matrix model introducing integration 
over an additional Hermitian matrix Y, required to have positive definite eigenvalues. We show 
how the Nambu-Goto version of the Green-Schwarz action is reproduced for the proposed model 
at the quantum level. The results are discussed in Sect. 4. Appendix A contains the proof of the 
N = 2 supersymmetry of the proposed matrix model for n — > oo. 



2 Interaction of branes in the IKKT matrix model 

In the large-n limit the matrices A„ and if) a become operators in a Hilbert space and the classical 
equations of motion ( |1.3| ) possess nontrivial solutions which possibly correspond to solitonic states 
in type IIB superstring theory. Among the solutions to eq. (|l.3| ), the distinguished role is played 
by the ones for which the field strength 

f,», = i[A ll ,A v ] (2.1) 

is proportional to the unit matrix. Only these classical configurations can preserve half of the 
super symmetries and thus can be interpreted as BPS states [l7L ||. This is why they are associated 
with D-branes of various dimensions. 

The solution which can be interpreted as D-brane of dimension p has the form 

Af = {B ,B 1 ,B 2 ,...,B p ,0,...,0), C' = 0, (2.2) 

where Bq, . . . ,B P are operators (infinite n x n matrices) with the commutator 

[B a ,B b ]=-ig ab l, (2.3) 

and a, b = 0, . . . ,p. Such solutions exist only for odd p, otherwise one can find a linear combination 
of B a 's, which commutes with all other operators. The solution corresponding to D-strings (p = 1) 
was studied in Q and was generalized to p = 3 and 5 in [R]] . 
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By a Lorentz transformation the skew-symmetric matrix g a ^ can be brought to the canonical 
form 

/ -u> x \ 
ui 



9ab 



V 



UJ p+l 

2 



2 





(2.4) 



For such g a b the operators B a form a set of I = (p + l)/2 pairs of canonical variables and (2.2) 
coincides with Ql.4|). In the coordinate representation they can be represented by 



Bn 



91, 



Bp-i 



iuidi, B p 



Qi- 



(2.5) 



The eigenvalues of the operators S a are uniformly distributed along the interval [— L a /2, L a /2] of 
the real axis, where L a /27r are compactification radii. In fact, we can assume that the support of 
eigenvalues covers the whole real axis, because L a should scale in the large n limit as n~i [?], ||. 
Each operator B a has n~ different eigenvalues, so the spacing between them, L a n~~, scales as 
rT~i . The product of the eigenvalue densities of the canonical conjugate variables £?2i-2 and Bii-i 
is fixed by the Fourier transformation, so we get 



2 

ni 



-^21-2-^21-1 



i 

n i 



(2.6) 



In this section we consider processes with interaction between p-branes of this type. We cal- 
culate first the interaction between two parallel p-branes moving with constant velocity, then the 
interaction between two p-branes rotated through some angle and finally the interaction between 
two anti-parallel p-branes (or brane-antibrane) . 



2.1 Scattering of parallel p-branes 

The multi-brane configurations correspond to the solutions of the equations of motion which have 
the block-diagonal form. Obviously, the matrices with two identical blocks describe a pair of 
superimposed p-branes. Let us first shift one of them by the distance b/2 and the other by —6/2 
along the {p + 2)-th axis. This results in the configuration of two parallel p-branes separated along 
the (p + 2)-th axis by the distance b from each other. Let us finally boost these stationary p-branes 
along (p+l)-th axis in opposite directions. Using the block-diagonal construction, the configuration 
of two parallel p-branes moving with constant velocity, with impact parameter b, is thus described 
by the following classical solution to eq. (|L1 

— 



A 

-tx, 



cl 



p+1 



A? 




a= l,...p, 



Bn sinh e / ' 



0, i=p + 3,...9. 



(2.7) 
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To simplify the calculation, we have chosen the frame where the two p-branes have opposite velocity 
v and —v along the (p + l)-th axis and 



tanh e . 



(2.8 



The interaction between the two p-branes to zeroth order in string coupling constant is deter- 
mined by the one- loop effective action of the matrix model in the background ( |2.7D , which have 
the form IM] 



W= ^Trln(P 2 V 



2iF, 



flV ) 



lTrln((^ + l^m (ii^l 



Trln(P 2 ) , 



(2.9) 



after the Wick rotation to the Euclidean space. The adjoint operators P M and act on the space 
of matrices and are defined by 



A 



cl 



F, 



(2.10) 



The block-diagonal form of the classical solution ( |2.7| ) shows that it is convenient to represent 
the matrices in adjoint representation as 2 x 2 matrices composed of re x n blocks. At infinite 
n these blocks become the operators acting in the same Hilbert space as B a . In the coordinate 
representation (2.5) they have the form of the functions of two sets of I variables — qj and g 2 - 
and B a act on them as derivative (left and right) and multiplication operators. From the definition 
of the adjoint operators, we find that P^'s act on Hermitian matrices as 



Pn 



Pi 



Pp-i 



Pn 



P 



p+l 



P 



p+2 



X Y 

yt z 

X Y 

yt z 



X Y 

yt z 

X Y 

yt z 

X Y 

yt z 

X Y 

yt z 



iui(dl +df)coshe ^ ^ ^ J , 



9?) 



X Y 

yt z 



iuji{d} + df) cosh e ^ * \ ^ 



(2.11) 



iui\ sinh e 



X Y 

yt z 



(d\ + d\)x (d\ - d\)Y 
-{d\-dl)Y^ -{d\ + dl)z 



o y 
-yt o 



(2.12) 



The only non-zero component of the field strength in the adjoint representation is Pi p +i which 
acts as 




2ui\ sinh e 




(2.13) 



This operator acts non-trivially only on Y and yt. The effective action fl2.9| ) vanishes for F^ v = Q, 
consequently we can trace only the action of P 2 on Y. The contributions of the eigenvalues, 
corresponding to the eigenfunctions with nonzero X and Z, are mutually cancelled. 
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Figure 1: The region of integration in the q + ,q -plane. 



Using the notation 



1 



! 1 to 1 t ,f) 



5f = -^(fl^db^), (2.14) 



we write the action of the operator P 2 as 

P 2 Y = j& 2 + 2J2 [(C) 2 " (^df) 2 ] - 2u\ cosh 2 e (d+) 2 + 2( 9 r) 2 - 2u, 2 sinh 2 e (<9f) 2 j F . 

(2.15) 

This expression is calculated in Euclidean space after a Wick rotation. We can clearly see that 
the last two terms that depend on q^ and give a Hamiltonian of harmonic oscillator, while 
all dj, and q~ for i = 2, . . . , I enter without their conjugate variables and can be simultaneously 
diagonalized with P 2 . Thus the eigenvalues of P 2 are given by 

1 1 
Eq,p,k = b 2 + 2 J2(qf + P* ) + 2 cosh 2 ep\ + 4wi sinh e (k + -) , (2.16) 

i=2 2 

where qi and p, are the eigenvalues of q^ and iuidf ', respectively. 

According to the discussion at the beginning of this section, the eigenvalues of the operators 
g*, g 2 are distributed from —Ln-xjl to Lii—\j2 with the constant density n l l l / L^i-i- This is 
not true for the eigenvalues of q~ , qf because the integration region changes under the change of 
variables (|2.14|) , as shown in fig. ||. As a result, the density of the eigenvalues qi in ( 2.16; ) decreases 



from y / 2n 1 ^/L2i_i at the origin to zero at rbl^i-i/V^- The scale of the variation is however of 
order L, which is negligible in the large-n limit. For convergent integrals the distribution of qi 
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and Pi can be taken to be uniform and equal to its value at the origin, y / 2n 1// '/L2i_i for and 
V / 2n 1 /7L 2i - 2 for Pi . 

Now we can utilize the results of Ref. ]?]] to bring the one-loop effective action for the given 
background to the form 



ajtl \ La J J k=0 



In 1 



16uj± sinh 2 e 
E q,p,k , 



1 ^-^ / 2siCc>i sinhe 

o E ln U- 



S1,...,S5=±1 
Sl...S 5 = l 



q,p,fc 



(2.17) 



The last term which originates from the integration over fermions can be rewritten as 



1 



E, 2sitJi sinh e \ , , / 4a;? sinh 2 e 



Sl,...,S 5 =±l 
S1...S5 = 1 



q,p,fc 



f?2 



(2.18) 



It is convenient to represent the logarithms in fl2.17 ) in the form of the integrals over a "proper 
time" s: 

oo 

u f ds 



ln - = / — [e 

v J s 
o 



( e -»*_ e ~ us ) 



The sum over k and the integrals over q and p then can be evaluated using the formula 



p 

TT \ 2 



£^ \2s/ 2 cosh e sinh (2wis sinh e) 



We finally obtain the following form 



i , _ i f°° ds f it\ 2 _ 6 2 s (cosh(4u;issmhe) — 4cosh(2cJissinhe) + 3) 



W = -n^ J] L- 1 / - - 



a^l 



« \S 



Defining 



coshe sinh(2(jJis sinh e) 
V p+1 = Y[L a 



p 

n 

a=0 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



and substituting (cf. eq. 



i 

n i 



n 



L2i~2L2i~l 
1-KU); 



^P+lII 



i=l 



27TL1J,- 



we now Wick rotate back to Minkowski space-time to obtain 



Vn 



1 f°° ds f 7r\2 _ 6 2 s (cos(4a;issinhe) — 4cos(2o;issinhe) + 3) 

coshe sin(2o;is sinhe) 



(2.23) 



(2.24) 
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where 

v 

V p = J] L a . (2.25) 

a=l 

It should also be noted that u>i ~ a' from dimensional analysis. We take 



2vrc/ , (2.26) 



for all i, which, as we shall now show, is the correct normalization to get agreement with super- 
gravity. 

We can compare the result we get from the IKKT model to that of Bachas [13] which is exact 
in b,v,a'. It is clear that the IKKT result does not agree with Bachas' calculation, for instance 
a comparison of the absorptive parts shows that they do not have the same poles (see below). 
There is, however, a regime in which the two results are identical. This is the regime in which 
supergravity, or alternatively the lightest closed string modes dominate the interactions. Thus this 
is a low-energy long distance approximation. This regime is characterized by b 2 3> a'. Bachas' 
expression for the one- loop effective action can be written in the more concise form [14|: 



poo 

F = VJ87r 2 a')- p/2 / dr r 
jo 



A+p/2 e -b 2 /2na'i 



OO 



-ttt/12 (1 — g- 2 ""™} 



n=l 



ie/n | it)) 4 



6>i(-i2e/7r | ir) 

(2.27) 

In the b 2 3> a' limit the integral is dominated by large values of r thus we can expand the theta 
functions in q = e~ nT to get: 

9\{z | it) = — 2g 1//4 sin(7r2i) + higher orders in q. (2.28) 

Keeping the lowest order in q yields 

F = -4iV p (8ir 2 a')- p/2 f° dr T -^ e -" 2 /^'r . (2 . 29) 



Py coshe ju 

By changing the integration variable to t = 1/r one can re- write the expression as 

,3 



4iV p (87r 2 a'r p/2 f°<& t 2 -P/ 2 e - fc2 '/ 2 ™' 

coshe Jo 



= -w r (^ r (2J " ,)3 ' P . (2.30) 

py ' coshe V 2 J b 6 ~P v ; 

One can compare this to the expression we have obtained above using the matrix model. If we 
perform a similar approximation in which b 2 3> a', we find that our expression is dominated by 
small s. We obtain an identical expression to Bachas in the above approximation scheme. Also, 
we mention that we get complete agreement between the matrix model and superstrings for small 
velocity and any 6^0. This is due to the complete cancellation of all the factors containing 
exp(— 7tt) in eq. ( |2.27| ). However, the next order term in the expansion of the small velocities does 
not agree. 

One might still ask how much the results from the matrix model deviates from the results 
obtained by Bachas. Due to the complicated integrals over theta functions, it is not so easy in 
general to give a quantitative estimate of these deviations. However, if we compare the imaginary 
part of the phase shiftEl, computed for D-branes by Bachas [13], with the same quantity computed 



2 The real and imaginary parts are defined by W — i(Re 8 + i Im 8). 
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from the matrix model, it is quite easy to see that in the limit where the velocity approaches light 



velocity, there is a very large physical difference between the results. The matrix result ( 2.24 ) has 



an imaginary part corresponding to poles coming from the trigonometric functions at the values 

Sk = 7, — -tt , for fc = 1,3,5, ... . (2.31) 
luii sinn e 

Computing the residue, we see (normalizing the parameters u>j as in ecL ( |2,26| )) that for small 
velocities we get a result which is the same as the one obtained by BachasH, 

Im s " WW < 2 - 32 > 

to leading order. It should be emphasized that this result is an independent check of the matrix 
model versus superstring calculations. This is because when we compute the large distance Re 6, 
the relevant region of integration is s — ► 0, whereas for Im 5 the relevant s is given by the position 
of the lowest pole. The result ( 2.32| ) also implies that the normalization ( |2.6| ) is correct for the 



imaginary part of the phase shift. 

However, for a large velocity the results differ. Introducing |1 



era ln(l-u)/2*J \n(s/M 2 p ) > 1, (2.33) 
where s stands for usual Mandelstam's variable and A4 P is the mass of the p-brane, we get 

/ _ \p/2-i 

Im S ra const [r^J exp(-b 2 M 2 p /s). (2.34) 

The main difference with the D-brane case is that these behave like black absorptive disks of 
logarithmically growing area, b 2 r ~ ln(s/A4p), whereas from the matrix model these black disks 
are much larger, corresponding to b 2 r ~ s/M. 2 . 

2.2 Rotated branes 

The configuration with two rotated p-branes can be obtained from the block-diagonal matrix with 
two identical blocks describing a pair of superimposed p-branes quite similarly to Subsect. fl2,l|). 
Shifting along the (p+2)-th axis by the distance b from each other and rotating in opposite directions 
in (p,p + 1) plane through the angle 6/2, we obtain the configuration of two rotated branes 



Ai = Ba y a = o,..., P -i, 

cl _ / B p cos f \ 

p { £ p cos§ J ' 

cl _ / 5 p sin| \ 

Ap+1 ~ { -Spsinf J' 

Af = 0, i=p + 3,...,9. (2.35) 



B a 



3 The units are a' = 1/2. Also, the velocity used here is one half of the velocity used by Bachas. 
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This looks just like an analytic continuation of eq. ( |2.7| ), Therefore all the formulas of this subsection 
are quite similar to those of the previous one. 

The interaction between these rotated p-branes to zeroth order in string coupling constant is de- 
termined by the one- loop effective action (2.9) in the background fl2.35|) . Repeating the calculation, 
we arrive at the Hamiltonian which has the spectrum (cf. ( |2. 16|) ) 

Eq, P ,k = b 2 + 2j2(q? + pf) + 2 cos 2 °- qf + Aun sin 2 °- (k + i) . (2.36) 



The final result for the interaction between two rotated p-branes, which are separated by the 
distance b we represent in the form 

2p 1 tt T 

w = _ 4n p+i — 7 n l- 1 



COS 2 



oo 



ds ( 7T 



x J — ( — J e b s tanh ^Wp+issin-J sinh^ ^Wp+issin-J . (2.37) 
o 

It can be obtained from ( |2.21| ) substituting e = i9/2. 

For large separation between branes we find, using the equality (2.6): 

/6-tA 4V v . -n- 1 sin 3 1 1 / 1 \ 

" r = - r (-5 £ )c5jr l/ *n^^fis=; + (p=;)- < 2 ' 38 > 

This expression correctly reproduces the supergravity result for the angular and distance depen- 
dence of the interaction energy between two rotated p-branes. An analogous formula for p = 1 is 
first obtained in @. 

2.3 Anti-parallel branes 

The interaction potential for two anti-parallel D-strings in the IKKT matrix model was calculated 
in 0. This calculation has been generalized to p-branes of arbitrary dimension |J. For complete- 
ness, we reproduce these results here using the same techniques as in the two previous subsections. 

The classical solution describing anti-parallel p-branes at the distance b from each other are 
represented by block-diagonal matrices 

= ( B Q a a = 0,...,p 

~ { -6/2 J' 

Af = 0, i=p + 2,...,9, (2.39) 
where B a and B' a obey the commutation relations 

[B a , B b ] = -ig ab l, [B' a , B' b ] = ig ab l. (2.40) 
The matrix g ab can be taken in the same form as in eq. Q2.4|). Thus we put 



B Q =iuj 1 d 1 , B 1 = q 1 , B p -. 1 =iuidi, B p = q h 

B'^iurft, B[ = - qi , B' v _ x = i^d u B' p = -q h (2.41) 
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where I = (p + l)/2. 

To calculate the one- loop effective action in the background ( |2.39 ), we perform the same steps as 
in Subsec. 2.1. First, we find all nonzero components of the adjoint operators defined by eq. (2.10): 



/ X Y 

°[Yi Z 

( X Y 

1 yt z 



Pp-i 



X Y 

yt z 

X Y 

yt z 



p A x Y 



X Y 

yt z 



( {q\-ql)X (ql + qf)Y 
{ -{<l\+<&)Y^ -(ql-Qf)Z 

iun(dt + $)[£ Y z \, 

( (qi-qf)X (qj + qf)Y 
{ -(ql+qf)Yi -(qj-qf)Z 

b ( -yt o ) 



(2.42) 



and 



-^2i-2 2i-l 



X Y 



-2ui 



o y 
-yt o 



yt z 

Again F^ u acts only on Y; the operator P 2 in this subspace has the form 

i 



i = l,...,l. 



P 2 Y = 6 2 + 2£[(g, 
I i=i 



+) 2 „ _ 



y . 



(2.43) 



(2.44) 



which coincides with the Hamiltonian of the /-dimensional harmonic oscillator. Therefore the 
eigenvalues of P 2 are labelled by the set of / positive integers and are given by 



£ k = 4X>i (h + - )+b 2 . 



1=1 



1 



(2.45) 



The form of the adjoint field strength ( |2.43j ) allows to use the representation for the effective action 
found in 0], the same as in eq. ( 2.17|) : 



W = nJ2 
k 



EM 1 



E 



In 1 



Si,— >*5=±1 
Sl...s 5 =l 



2£ UiSi 

i 



(2.46) 



Neither the quantities q i , nor their conjugate variables enter ( 2.44] ) , so the trace over them gives 
an overall factor of n. 



Using the same proper time representation for logarithms as in eq. (2.19) and the equality 

„-b 2 s 



\\ 2 sinh 2u>iS ' 



(2.47) 
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we obtain for the effective action: 



W 



-In 



ds 



o 



y^fcosh 4u>jS — lj — A^Y[ cosh 2ujiS — 1 



n 



i 



. 2sinh2o;,-s 



(2.48) 



For large separation between the branes, this result reduces to 

/ \ 2" 



w 



1 

16 



n r 



7 — p 



UJ.; 



n 



i /2\ 7 ~ p 



b 



+ 



b 9-p 



(2.49) 



Equations ( p. 48 ) and ( |2.49 ) coincide with the previous results |7|, ||] obtained by a slightly different 
technique. 

These results of the matrix models are to be compared to the superstring calculations which are 
given in the open-string language by the annulus diagram. The superstring result for the interaction 
between anti-parallel Dp-branes reads jy], 14] 



W 



-V, 



P+i 



dt 



1 



p -b 2 t/2na'-l n^=l(l - Q 



2n-l\8 



(2.50) 



o ^(SvrVt)^ riT=i(l-^) 8 

with q = e~ nt . We see that the superstrings and matrix-model answers agree only at large 
distances^ quite similarly to the cases of moving and rotated branes. This suggests to modify the 
matrix model to better reproduce the superstring calculation. 



3 The NBI-type matrix model of IIB superstring 

In the IKKT model the matrix size n is considered as a dynamical variable. The partition function 
includes for this reason a summation over n: 

oo . 

Z=J2 VAVil) e~ s , (3.1) 

71=1 ^ 

where the action is given by eq. (|1.1|). This construction is proposed in ]7| as the matrix-model 
analog of the Schild formulation of type IIB superstring given by the path integral 

Z = J V^jVXVip e" 5schild (3.2) 



with the action (1.2). 

In this section we propose a modification of the IKKT matrix model which appears to be a 
more analogous to eq. (^), and which reproduces the Nambu-Goto version of the Green-Schwarz 
superstring action after the integration over the introduced additional Hermitian matrix Y %3 with 
positive definite eigenvalues, which is analogous to ^fg in eq. (3.2). The classical action has the 



form (|l.5| ), which yields the following classical equation of motion for the y-field: 

- (y-^Afy- 1 ).. + p8ij = o, (3.3) 



4 Tseytlin jl^] has conjectured an alternative interpretation of the classical solutions in the IKKT matrix model 
as D-branes with magnetic field, in analogy with previous work |i| on the matrix model J5J. We do not discuss such 
a point of view in this paper. 
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whose solution reads 



Y = - ] I^-[A„A^. (3.4) 

Here — LA^,^] 2 is positive definite, since the commutator is anti-hermitian. The square root in 
( |3.4| ) is unique, provided Y is positive definite which is our case. After the substitution fl3.4|) , the 
classical action ( |1.5| ) reduces to the classical action of the NBI (non-abelian Born-Infeld) type ( |1.6[ ). 

In this section we show that even at the quantum level, it is possible to modify the classical 
action ( |1.5| ) such that we obtain the Nambu-Goto version of the Green-Schwarz type IIB superstring 
action. 

3.1 The Yang-Mills dielectric matrix model 

Let us consider the matrix model defined by the action 

Se = ~\ Tr (^LV A,] 2 ) + V(Y) " f Tr {^[A^]) , (3.5) 

where Y is a new n x n Hermitian matrix field taken to be positive definite, and V(Y) is a 
"potential". For reasons which become clear later, the potential is taken to be 

V(Y) = (3 Tr Y + 7 Tr In Y. (3.6) 

The partition function is then given by the functional integrals! 

Z e = J VA^ £>V> VY e~ S *. (3.7) 

As discussed in the appendix the action is invariant under 

= l{Y-\[A„A u ]}T^e, 

<5 (2) V> = f, 

S^Ap = (3.8) 

in the limit n — > 00. The field Y is assumed to be invariant with respect to this transformation. 

The action ( |3.5| ) differs from its classical counterpart ( |1.5| ) by the second term on the right-hand 
side of eq. (|3.6j ). We can e.g. associate this term with the measure for integration over Y rather 
than with the classical action. The classical action ( |1.5[ ) can be obtained from Q3.5| ) in the limit 
a ~ f3 ~ 00, a/ [3 ~ 1. This limit is associated with vanishing string coupling constant since 
a ~ (3 ~ g~ l , i.e. with the usual classical limit in string theory. The matrix model with the action 
( |3.5| ) can be considered as the large n reduced model for a ten-dimensional non-abelian "dielectric" 
theory of the type introduced by 't Hooft |l7j several years ago for the abelian case. In this picture 
the quantity 1/Y is the dielectric function e(Y), which is then governed by the potential V(Y). 
Although the present matrix model and its non-reduced counterpart are much more complicated 
than the abelian version!, we shall call the model given by eq. ( |3.5D the Yang- Mills dielectric matrix 



5 In order to have a non-trivial saddle point for n — > 00, we need to assume that the positive constants a and (3 
are of order n. 

6 This applies to the physics as well as to the mathematics. 
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model. This is the reason for the notation S e for the action. Alternatively, one could interpret the 
field Y as a rather rudimentary metric. 

We start by doing the Y — integral. Since the tp— dependent term in eq. ( |3.5| ) is independent of 
Y, it is sufficient to consider the integral 



F{z) = J VY exp (-j Tr ~ P Tr ^ - 7 Tr InY 



(3.9) 



Here z = —[Ap,Ai,] . The integration over the "angular" variables in eq. (|3.9|) is of the Itzykson- 



Zuber type |18[ , and we therefore get 
F{z) = n\ a^- 1 ^ 2 TT pi TT / d Vi 



A 2 (y) 



A(l/y)A(z 2 



exp -a zf/Ayi - (5 ^ y { - 7 ^ lny^ 



where A is the Vandermonde determinant 



A(x) = TT(xj — xj) = det x 

ki 

i>3 



k-1 



(3.10) 
(3.11) 



In eq. (|3.10|) the quantities zf and yi are the eigenvalues of —[A^^Ay] and Y, respectively. The 
unitary matrix which diagonalizes Y has thus been integrated over as in the Itzykson-Zuber paper 

Hi- 

Now eq. ( [3.10 ) can be written 



n— 1 n 



A(z 2 )^(z) = n! a^" 1 )/ 2 J] ^ II / dVi V?' 1 Ufa ~ Vi) 

p=l i=l J ° i>j 

x exp I -a z i/ 4 Vi ~ $ Vi ~ 7 J2 ) ■ ( 3 - 12 ) 

We shall now choose 7 in such a way that the result of the yi— integrations gives a result which is 
"as string-like as is possible". As we shall soon see, this amounts to taking 



7 = n - 1/2. 



(3.13) 



Using the well known Bessel integral (from an integral representation of i£_i/2 and the explicit 
formula for this function) 



dy 



y 



1/2 



-az 2 /4y-Py 



'a/3 z 



we obtain by use of ( 3.13| ) 



n-l 

n! a"( n - 1 )/ 2 TT p\ det(- 
1 ki 



\Jtt/P exp {—\[afiz. l 



The matrix whose determinant is to be calculated reads 

d 



A ki = (-1 



,fc-i 



d(3 k 



(3.14) 



(3.15) 



(3.16) 
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or, explicitly, 



A 2i 



A 



1 ,— 
1 



I afizi 



I a(3zi 



+ 



1 

2$ V/3 



'apzi 



3/ 



P 2 - 



lafizi 



+ 



4 (3 2 



-Zi e 



TT 



+ 4/3 2 V /3 6 



(3.17) 



The second term in the expression for An is proportional to the first line of the matrix, Au, and 
can be omitted in the determinant. The same property holds for all lines of A - only the result 
of the differentiation of the exponential survives in the determinant, while the terms coming from 
differentiations of the various pre-exponential factors are linear combinations of the previous lines 
of the matrix. These factors include e.g. the term in A^i which is linear in Z{. This term is 
proportional to the first term in An, etc. etc. 
Hence, 



A(z 2 )Hz) 



C det 

ki 



fc-1 



C AO) exp 



(3.18) 



where the constants C and C are given by 



n-l 



n-l 



C = n\ a n{n ~ 1)/2 J] P! and C = n\ ((a) 3/2 /2 V / ^) n(n_1)/2 (^//3) n/2 II P l ■ ( 3 - 19 ) 
p=i p=i 

In these equations, always means the positive square root of the positive quantity zf. 

Eq. ( |3.18| ) is one of the main results of this subsection. It shows that after the exact integration 
over Y, for any n, the action becomes linear in the variables Zi, although the original action S e in 



eq. (3.5) is quadratic in these variables. To the best of our knowledge, this result is a new matrix 
model result, and it may therefore be of interest also outside the present framework. Expressed in 
terms of string language, we shall find that we have the option of obtaining the (supersymmetric) 
Nambu-Goto type of action, as we shall discuss shortly. The result ( |3. 18|) is possible because of 
the choice ( 3.13| ) of the power 7, which allows us to use the well known explicit (exponential) form 
for the Bessel function K_i/ 2 - However, it should be noticed that we still have the ratio of the 
Vandermonde determinants in eq. ( p.!5| ), 



A(z)/A(z 2 ) = l/II( 



(3.20) 



i>j 



To proceed from eq. ( |3.15| ) we now show that we have the identification 



Tr 



'[An, A v ] 



(3.21) 
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where it should be remembered that the zf's are the positive eigenvalues of —[A^,A^. Following 
Dirac we define the square root of a matrix by its formal power series. Thus we write 



i p=0,l<p 



P 
I 



p-ljl 

i > 



(3.22) 



where c p = 1 • 3 • 5 • ... • (2p — 3)/2 • 4 • 6 • . . . • 2p. Now let U be the unitary matrix which diagonalizes 
the square of the commutator, 



diag(>?,2f 



! • • • i "n ) 



-XJ\A^A V \ 2 XJ\ 



(3.23) 



Then we can write in eq. (|3.22j) 



Q = )("!)% Tr (u[A (tl ,A Vl ]Wu...U[An,A Vl ] 2 rf 

= Tr£c p (-[A„A v ] 2 -l) P = TryJl + i-lA^Arf-l) 



TrJ-[A„,A v } 2 . 



(3.24) 



This verifies the identification ( 3.21[ ). 

By means of the result displayed in eq. ( |3.24 ), we can rewrite the partition function (3^) by 
use of eqs. (EOD and (|1D in the following way 



Z £ = J VAfjV^VY exp (| Tr (± [A^ A v ) 2 ^ - (3 Tr Y 



I a 
(n--)Trlny + -Tr (^[A^]) 



a 



C / [VAMexp [-^TvJ-{A^A v f + -TT{^[A^\) 



where the measure in the last expression is defined by 

VA^ 



\VA^\ 



(3.25) 



(3.26) 



The square root of a matrix can alternatively be defined by the integral representation 

c(o+) 



V™ J-, 



dtr 3/2 e- tlA - A » ] 



where the contour of integration encircles counterclockwise the negative real axis which provides convergence. Thus 
we write 

r(0+) 



which can be represented as 



(0+) 



dtr 3/2 Tve- t[A ^ Av? =Tr V / ^[I^XF 



All the formulas here and below can be rigorously derived using this representation. 
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The result ( |3.25| ) shows that the "dielectric" action S e is equivalent to an action which can be 
considered as a strong coupling non-abelian Born-Infeld modefl The new action in ( |3.25 ) is no 
longer invariant under the transformations (|3.8|). The first of these should be replaced by 

S^^ = ^{(-[A a ,A p ] 2 r 1 / 2 ,[A„A v }}r^e. (3.27) 
The other transformations in eq. fl3.8|) are unchanged. Thus the effect of performing the Y- 



integration is not only to produce a new action, but it also produces a new measure ( 3.26 ) and a 
new transformation property ( |3.27|). 

The anti-commutator in eq. ( |3.27[) is superfluous, since any operator commutes with its own 
square. However, it follows from the appendix that the action in the second equation ( 3.25| ) is only 



invariant under the transformation ( 3.27| ) in the limit n — » oo. This is because the last two terms 



on the right hand side of eq. ( |A.4 ) do not vanish when Y 1 is replaced by the inverse square root, 
as in ( p7D - 

3.2 Connection between the dielectric matrix model and superstrings in the 
large n limit 

In this section we shall discuss the connection between the super- "dielectric" matrix model and 
superstrings. The possible connection between the large n limit and strings have been discussed 
by many authors. We shall use in particular the QCD-discussion by Bars in ref. piJl , where further 
references can be foundtl The starting point is an expansion of the ^4^-field in terms of SU(n)- 
generators /k,k = (ki,k2), 

(^' = ^E fl J(Wi'' ( 3 - 28 ) 

k 

where are expansion coefficients to be integrated over in functional integrals, and where C\ 
is some normalization constant. The expansion ( |3.28| ) is to be compared to the corresponding 
expansion of the string variable X^(a,r), where a and r are the usual world sheet variables. In 
this case we have 

X>,r) = C 2 ^aJe iW , (3.29) 
k 



where a = (a, r). The Z-matrices can be constructed explicitly in terms of Weyl matrices [21|, and 
they satisfy the commutation relation 

Tl 

fciA 2 ] = sin(2vrki x k 2 /n) q< 1+ k 2 -» «( k i x k 2) Zki+k 2 , for n -> oo, (3.30) 

where ki x k 2 = (^1)1(^2)2 — (^2)1(^1)2- I n the following we shall always perform the n — > 00 
limit in the way done in eq. ( 3.30| ), even inside sums over the k's. We do not know a rigorous 



justification of this. Obviously there is a good chance for the validity of this approach if infinitely 
high modes (n — > 00) are not dynamically relevant. For the bosonic string this is probably not 
true, since this string oscillates infinitely much at short distances due to the tachyon instability. 
However, the assumption may be correct for stable superstrings. 



8 The reader should recall that the abelian Born-Infeld action |19| has the Lagrangian £ = -cy/l + F^/2c, where 
c is a constant. In the strong field limit C behaves like —yj c/2^J ~F%,. It was suggested many years ago that the 
strong field limit of the abelian Born-Infeld should give a field-theoretic description of strings J2^] . 



We mention some of these works in 
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The results for SU(n) mentioned above can be compared on the torus with the area preserving 
diffeomorphism. This is discussed in the paper by Bars [21] (for a discussion of SU(n) on the sphere, 
see Floratos, Iliopoulos, and Tiktopoulos [23]), and we shall not repeat this discussion. We only 
wish to mention that this approach corresponds to taking into account the local subalgebra, but 
ignoring the global translation generators. The central extensions are thus ignored. 

Let us again consider the quantity Q defined in eq. ( 3.22| ). Inserting the expansion ( |3.28| ), we 
get from eq. (|3.24j) 



Q 



Tr E 

P ,l 



P 
I 



(-1)% 



-Cf a " a ™ a £°£( n x m )( r x s ) ln+mh 



E 

p,i 



p 
i 



nmrs 

/ 



Cf J2 ]J (" >< m)(r x s) 



nmrs i=i 



(3.31) 



3=1 



In the limit n — > oo the trace of a product of Z's with different indices produce a Kronecker delta§ 
(up to a normalization factor). Therefore eq. ( |3.31 ) becomes 



Q 



p 
i 



n E 

P ,i 

x "ni+mi+ri+si...+n i +m i +r i +s i ,0 



21 I 
nmrs j=i 



ir*<* (^cf) E n(^«<:( n x m )( r >< s 



(3.32) 



It should be emphasized that this simple result is valid only because the n — » oo limit is taken. 
Otherwise the trace in eq. ( |3.31| ) yields a more complicated result than displayed above. This 
expression can be rewritten by use of the expansion ( 3.29| ) for the string coordinates, 



Q 



n 



2lT 



(2vr)2 
n 

(2^)2 J u 
1 

2(2^)3 V C 2 



jf^'E ( \ J (-l)^((n/47r) a (C a /C 2 ) 4 {X»,X»} 2 PB 



dV^/l + ((n/47r) 2 (Ci/C7 2 )4 {X"(<r),^(a)} 



(3.33) 



Here {a, 6}pb is the usual Poisson bracket of a and b. 
Inserting this result in eq. ( |3.25| ) we finally obtain 



Z e ^C [VX^] exp 



la 



d l a ( Vaf3j{X^a),X u (a)} 2 PB - ^-e ab d a X^F^ 



(3.34) 



This follows from 



Tr l m lr = (n 3 /(47r) 2 ) S m+r n 



and repeated applications of the relation lmlr = (n/Air) exp(27ri(m x r)/n) lm+r- The exponential factor produces 
one plus terms of higher order in 1/n. Thus, 

Tr Zm 1 .../m s = n(n/47r) s <5m 1 +...+m s ,o 



to leading order. 
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The last term in the action follows from expanding ip in a form similar to (3.28) and using that for 
n — > oo 



Tr^rq^,^] = £ a^ r r^ m (n x m) Ttl r l n l m 

d 2 ae ab d a X^TM, 



8(2vr) 5 C 2 



(3.35) 



where is the relative normalization of the ^-fields. The coefficients in ( [3.33|) and ( |3.35| ) have 
been absorbed in a redefinition of a and f3 in the result ( |3.34 ). Thus, e.g. 



(nCi/C 2 ) 2 /2(2^) 3 v / o^ -> yfcp. 



(3.36) 



It should be emphasized that the measure in the result ( |3.34| ) is defined through (3.26). Also, the 
transformation property ( 3.27 ) must be replaced by 



5 dty = 1 

2 y{ x a,X /3 } 2 PB 

in order to ensure invariance of the action in (3.34). 



{X^X u } PB Y^e 



(3.37) 



4 Discussion of the results 



The main result from the dielectric matrix model is given by eq. ( 3.34 ). We shall now discuss this 
result. First, one might wonder if this formula cannot be applied with ip = 0, so that the bosonic 
string would emerge from the dielectric model with ■0 = 0. As already discussed in connection with 
eq. ( 3.30 ) this is highly unlikely, since we perform the limit n — > oo inside sums, like e.g. in the 
transition from (3.31) to ( 3.32j ). This is allowed if the infinitely high modes are not dynamically 
relevant. However, we know that due to the tachyon, at a finite distance of the order the square root 
of the string tension, the bosonic string becomes unstable, due to the relevance of infinitely high 
modes. These causes the area of the world sheet to become infinite, due to an infinitely oscillating 
string. Hence it is not permitted to interchange the sum over modes and the limit n —* oo. For 
superstrings, the situation is much more hopeful, since it is stable without a tachyon. So although 
we do not have a mathematical proof that the summation over modes can be interchanged with the 
limit n goes to infinity, there are physical reasons to believe that this is possible for superstrings. 

The comparison of the results obtained from the matrix model with the ones from superstring 
theory provides a check to what extent the superstring theory can be described by the matrix 
model. The calculation of the interaction between D-branes is one of such checks. It is natural 
to think that such calculations in the matrix model correspond to loop expansion around certain 
large- n classical solutions. The one-loop calculations of Sect. 2 are performed without summation 
over n, as proposed in [f7|, or without integration over Y, as proposed in the previous section. 

It is worth mentioning, however, that the classical solutions (|2.2j ), which are associated with 
D-brane configurations, are also classical solutions to the NBI-type matrix model. The reason is 
that these classical solutions are BPS states and the commutator [A^, A u ] is proportional for them 



to the unit matrix. The same is true for the classical value of Y, as it follows from eq. (3.4), so the 
classical equations of motion of the NBI model: 



A^{Y-\[A^A U ]}] =0, [A^(T^) a ]=0, 



(4.1) 
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are also satisfied. 

A more general property holds in the large-n limit when any classical solution of the IKKT 
model is simultaneously a solution of the classical equations of motion of the NBI model. To show 
this, let us rewrite the equations of motion ( |4.l| ) and ( |3.4[) for bosonic matrices in the form 

{[A» ,¥-*], [A„A U ]} + {Y-\ [A», [A^A v ]]} = 0, Y 2 = [A„A U ] 2 . (4.2) 

For a solution of the IKKT model, the second term on the left hand side of the first equation equals 
zero. At infinite n, when the commutators can be replaced by the Poisson brackets, the first term 
also vanishes, since the large-n classical equations of motion imply d a Y 2 = = d T Y 2 |24|] and thus 
the Poisson bracket {A im ,Y~ 1 }pb is equal to zero. 

In this paper we have not discussed the large n saddle point configuration of the integral over 
the matrix field Y in the partition function Z e . It should, however, be emphasized that such a 
calculation is very different from the corresponding "classical" saddle point calculation, valid for 
a ~ j3 — > oo. In the large-n saddle point, the logarithm of the Vandermonde determinant enters, 
and one needs to determine the spectral density of the eigenvalues yi which in turn determines the 
value of the commutator [An , A v ] . This is most easily seen by summarizing our result in the form 

Z e = J VA^VY J[( Zi + Zj ) exp (~ Tr A u ] 2 ^j - TrY 

_( n _ I)Trlny+|Tr (^[A^})^j 

= Cj VX^Vi, exp - J d 2 a U^yf{X^), X u (a)} 2 PB - % -^e ab d a X^T ^d b ^j ,(4.3) 

where now the measure in the last functional integral is the standard one, whereas the corresponding 
quantity in the first functional integral is not, due to the factor Y\.( z i + z k)- This implies that in 
evaluating the large n saddle point this additional z dependent factor in the measure should be 
taken into account. Equation ([O]) for the A^-field with a nontrivial distribution of the eigenvalues 
of Y possesses undoubtedly a richer structure than eq. (|Q|). 



The matrix model given in eqs. (3.5)-Q3.7|) can presumably be considered as a large n reduced 



Eguchi-Kawai model for the field theory with the action 

S fie i d oc J d 10 x Tr (jY^F 2 , + ^VDrf + V(Y)j , (4.4) 

where V(Y) is the non-polynomial potential given in eq. (|3.6|). As usual, we have 

Ffu, = dpA v - d v A^ + i[Ap,A v \, and = + (4.5) 

Because of supersymmetry for n — > oo, we do not expect that quenching is necessary |25|, [?J. 
Therefore the non-polynomial action (|4.4j ), in the limit where n approaches infinity, could be 
considered as an effective field theory for superstrings. 

From this point of view, one could consider the field theory as a regulator for the Green-Schwarz 
superstring. For example, this might be useful in numerical simulations. However, one could also 
ask if the underlying field theory could be of direct physical interest. A possible scenario could be 
the following: Suppose that at the Planck scale or below there exists a description in terms of some 
unified field theory (probably non-polynomial) with a high order group, like e.g. SU(n) with n very 
large. Such a theory could then, as exemplified by our discussion above, effectively be equivalent 
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to a Green-Schwarz superstring theory. There would therefore exist a dual description of the very 
early universe, either as some unified field theory, or as a string theory. Of course, it goes without 
saying that many problems should be solved, before such a dramatic scenario can be said to be on 
a satisfactory basis. 
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Appendix A On the supersymmetry of the Y-integral 

In this appendix we shall show that the action (|3.5|) is invariant under the symmetry transformations 



( [3.Sj ) in the limit n — > oo. If we apply the transformation (|3.8| ) to the action (3.5), we obtain after 
some calculations 

\ Tr (l[A^A v \^ - I Tr (e m (T,T ,) mn ^ n [A^{[A^ A^Y- 1 })) , (A.l) 



as well as 



Tr(^L4^]) - \ Tr (^T 6 )^ ^.{^^.r 1 }]] 



1 _ „ % 

2 

* Aiafi\\...\-i Tr U (pOpllpA* r A 7 ) £p [^ ){[ ^ Q) A^]^- 1 }]) . (A.2) 



7! 4 



Here the quantity {a, 6} denotes the anti-commutator of a and 6, and should not be confused with 
the Poisson bracket. In deriving this result we used the expansion 

p/i-pa/3 _ ^£tap/3 _ ^M/3p« ±_ iMxp\i...\r pllpAi pA 7 ^ g-j 



If Y was a c-number, the last term in ( |A.2| ) would vanish for symmetry reasons, and the first term 
in this equation always cancel the expression on the right hand side of (A.l), corresponding to the 
well known invariance of supersymmetric Yang-Mills matrix theory. However, the presence of the 
non-commuting Y makes life harder. Here we shall show that the last term in eq. ( |A.2j ) vanishes 
in the limit n — > oo. Using the expansion 

^{{A^AplY- 1 }] = [A^A^ApWY' 1 + Y- x [A ll ,[A ol ,Ap]\ 
+[A a ,A /3 ][A tl ,Y- 1 ] + [A tl ,Y- 1 ][A a ,Ap], (A.4) 



the first two terms give zero contribution when inserted in the last term in eq. ( |A.2| ). The critical 
terms are thus the last two terms on the right hand side of (A.4). Consider one of these terms in 
the large n limit, 

mnpr 



[A a ,Ap][A^Y x ] = J2 ^ r ^ sin (~ m X n ) 2^ sin ("^" P X r ) /m + n W> ( A - 5 ) 
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where we used the expansion ( 3.28[ ) and the commutator ( p.32| ) , as well as an expansion of Y 1 , 

V 1 )) = E(y _1 ) r (h)i- ( A - 6 ) 

r 

Taking the limit n — > oo and using the expansion of ifj we get 

Tr(V [A^ApU^Y- 1 ])^ J2 ^ S C «" (y _1 ) r (m x n) (p x r) Tr(y m+n / P +r) 

smnpr 

-»■ £ V' 8 a™ «" «m (y _1 ) r ( m x n ) (P x r ) 4+m+n+p+r,0 
smnpr 

= J ' d 2 a^ {X a (a),Xp(a)} PB l/\f^a)} P B, (A. 7) 

using the expansion ( [3.29D . Also, we used 



VVff(^) = EGr 1 ) 1, e?r<T - ( A - 8 ) 

r 

Now we have 

e" " 3 - {A^, l/VsWA a , A>} PS = 2e^»- (j^A^l/Vs)' - A^A^l/Vs)') , (A.9) 



where dot and prime denotes derivatives with respect to r and a, respectively. The expression 
on the right hand side of ( A.9 ) is easily seen to vanish, since the first term inside the bracket is 
symmetric in fi and a, whereas the last term is symmetric in fj, and (3. Thus, for n — > oo we have 

Tr^iA^A^A^Y- 1 ]) - e^- J d 2 a i> {A /t , l/y/g} PB {X a , X p } PB = 0. (A.10) 



The proof of the cancellation of the cubic ip terms, which emerge in the action (|3,5| ) under the 
symmetry transformations ( |3.8| ), is the standard one. 
Hence the action (|3.5| ) is supersymmetric for n — > oo. 
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